On 



Lattice Approximation of Quantum Statistical Traces at a 

Complex Temperature 

J. Lukkarinen a ) 

Helsinki Institute of Physics, P.O.Box 9, 00014 University of Helsinki, Finland 
We prove that the simple condition on the potential V, J exp(— tV) < 00 
for all t > 0, is sufficient for the lattice approximation of Tr (Ae _/3 ^) with 
Re [5 > to work for all bounded functions A and a large class of potentials. 
As a by-product we obtain an explicit bound for the real-temperature lattice 
kernels. 
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O ■ I- INTRODUCTION 

in 

Ever since the introduction of the path-integrals by R. Feynman, the path-integral for- 
mulas with a complex "temperature" have presented an elusive mathematical problem — for 
instance, the Wiener measure, which is so useful for positive temperature statistics, cannot 
be extended to the complex case.0 In fact, since real time processes and real temperature 
statistics can be identified with each other via analytic continuation, the present numeri- 
cal non-perturbative methods in quantum field theories rely mostly on computations with 
t^. I lattice approximations to the imaginary time, i.e. real temperature, path-space measures. 

There have already been proposals for numerical evaluation of the complex path- 
integrals@i, but only now can we establish rigorously the assumed relations between the 



complex lattice integrals and the continuum traces. Also, the result of Theorem III.5 has 



been useda as an assumption in the derivation of a lattice approximation to the quantum 
microcanonical ensemble and the theorem now yields a wide range of potentials for which 
this approximation is valid. 
^ ' For an introduction to the mathematics of operators on a Hilbert space see e.g. the 

second part of the comprehensive series! by Reed and Simon. The path-space measures and 
their applications are discussed in detail in another book! by Simon and the proofs of the 
statements made in the following section can be found from Section 6 of that book. 

II. NOTATIONS AND DEFINITIONS 

We will now consider quantum mechanics on the separable Hilbert space L 2 (M. n ); in the 
following, n will always denote the number of dimensions of the parameter space and we 
will use boldface letters for vectors of M n . For simplicity, we have also adopted the following 
non-standard notations: -Dx denotes the n-dimensional Lebesque measure with respect to 
the variable x and D N x stands for the product measure Ylk=i -^ x fc- Operators will be 
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distinguished by the "hat" — especially the multiplication operator defined by a measurable 
function F will be denoted by F. 

The free Hamiltonian is defined in the standard way, 

T) 2 ft 2 

H ■= — = ——V 2 

2m 2m ' 

with the usual domain for this operator. To simplify the following discussion we will from 
now on use the natural units h — 1 and, in addition, consider only the case m = 1. The 
results for the situation m^l are easily reproduced from the following formulas by simply 
replacing t by t/m and the potential V by mV . 

It is well-known that Hq generates an analytic semigroup on the right half-plane, where 
e~" H ° are bounded integral operators for Re (3 > and they are given by the integral kernels 

F(a ' M):= ^ exp( -^ |a - b|2) - 

It is also evident that these integral kernels satisfy the semigroup property exactly, i.e. that 
for all a, b G K" and Re (3 U Re f3 2 > 0, 

Dc P(a, c; ft)P(c, b; f3 2 ) = P(a, b; fa + ft). 

The potential V is assumed to be locally in L 1 (M n ) and positive, but the results easily 
generalize for all real L\ oc potentials bounded from below (just replace V by K+ess sup(— V)). 
The Hamiltonian H = H + V can now be defined as a sum of quadratic forms, when it is 
self-adjoint and bounded from below and, therefore, e~^ H (as well as e~^ v ) is an analytic 
semigroup on the right half-plane. In addition, we shall also require that 

jDae- tv(a) < oo, for alU > 0, (1) 

so that e - *^ is trace-class for t > by the Golden, Thompson, Symanzik -inequalityB. 

Under these assumptions e~ tH is an integral operator for t > and it has an integral 
kernel K^, which can be obtained with aid of the conditional Wiener measure fi: for all a, 
b G R" and for alU > define 



In addition to ([T]) we shall now require the potential V to be such that the Wiener integrals 
can almost everywhere be obtained from the lattice kernels 



Nn N . N 



K N (z,b;t) := [d n -^(^-] * exp[ - |£ |x fc _ x - x fc | 2 - 1 £ V(x k ) 



x =a ' 
J xjv=b 



i.e. that 



^^(a, b; t) = lim K N (a, b; t), for alH > and for almost all a, b e R n . (2) 

N^oo 
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This requirement is clearly satisfied for all a and b if V is continuous, since then jdsV{uo(s)) 
can be approximated by suitable Riemann sums leading to @ by using the definition of the 
Wiener measure. 

We shall also need the operator 

f{{3) := e-P&e-P 9 , 

which is clearly a bounded integral operator for Re/3 > 0, since then ||T(/?)|| < 
llg-^ojjiig-pVjj < \ Notice that Kn(-,-;/3) is nothing but the integral kernel of the op- 



erator T{jA . 



III. TWO LATTICE THEOREMS 



In this section we will prove two theorems concerning the convergence of traces of lattice 
operators. The first theorem will show that also certain unbounded multiplication opera- 
tors are measurable using the (real-temperature) Gibbs statistics, if the potential increases 
sufficiently fast at infinity. This will be proven by using the explicit bounds given by the 
following lemma and its two corollaries. The rather technical proofs of the lemma and the 
corollaries are included in Section |V| — the lemma itself is an extension of a result derived 
by SymanzikP for the continuum kernels. 



Lemma III.l For all a, b in ~&L n , for all t > and for all positive integers N, 

K N {a,b;t) < C n t^ jDcP( a ,c;t)e- tv ^P(c,b;t) + j^^pi ) 

+ _L e -*nb)p( a)b;t)) 

where C n depends only on the dimension n. 



(3) 



Taking the limit N — > oo in the previous lemma shows that (W is valid for almost all a, 
h and for allt > also when N = oo. 

Corollary III. 2 Assume that in addition to ([Jj the potential satisfies 

^a|a| r e-* y(a) < oo, 



for some r > and for all t > 0. Then there exists functions F r and G r , continuous on 
(0, oo), such that for all t > and for all N < oo, 

J Da \a\ r K N (a, a; t) < F r (t) 

and for all N , including oo, 

jDaDb\a\ r K N (a,b;t) < G r (t). 
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If the above assumption on the potential is true for and r then it is clearly true for all 
values s in (0,r) by the inequality |a| s < 1 + |a| r . If the potential has no singularity in a 
neighborhood of the origin, it is enough to find one r > for the statement to be true for 
all s with < s < r. 

Corollary III. 3 For all N, including oo, and for all t > 0, ess sup K^(a, b; t) < C' n t~% . 

a,bGM™ 

The convergence properties of the lattice traces are easily derived from these results; we 
have only to be careful that the limit of the traces of the integral kernels does give the trace 
of the continuum integral operator. This leads to the appearance of the factor of two in Eq. 
(|j), which will be discussed more after the proof of the theorem. 

Theorem III. 4 Let A be measurable function bounded by a degree r > polynomial and 
assume that the potential additionally satisfies 

y_Da|a| 2r e-* y(a) < oo, for all t > 0. 

Then for allt>0 and N > 2 the operators Ae~ 1 ^ and AT(t/N) N are both Hilbert- Schmidt 
and trace-class and 

T r (i-)=ltaT r (if(±)"), (4) 



Tr [AT{t/N)") = J .Da A(a)ifjv(a, a; t). (5) 
Proof: The requirement for A is equivalent to finding a C > such that |A(a)| < C(l + |a| r ) 



for all a. It follows now from Corollaries III. 2 and III. 3 that for 1 < N < oo 



y > L)aDb|A(a)K i v(a,b;t)| 2 < oo. 



The operators Ae tH and AT(t/N) N are given by the integral kernels v4(a)iCv( a ; b; t) and, 
therefore, they are Hilbert-Schmidt. 

Since Ae~ tH is a product of the two Hilbert-Schmidt operators Ae~^ H and e~^ H , it is 
trace-class. Moreover, the trace can be expressed in terms of the integral kernels of these 
Hilbert-Schmidt operators: 

Tr(Ae-^ = jDaDb A(a)lf 0O (a, b; ^^(b, a; *-). 



Using Lemma |111.1| and Corollary 111.3 to give an upper bound for the Lebesque dominated 



convergence theorem shows now that 

DaDbA(a)if (X> (a,b;-)iif 00 (b J a;-) = lim DaDb A(a)K N (&, b; -)K N (b, a; -] 

2 2 N—>oo I 2 2 



lim Da.A(a)K 2 N(a.,a.;t), 

N^oo J 
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where the second equality follows from Fubini's theorem and a careful inspection of the 
definition of K 2 n- Therefore, 

Trfle"* 5 ) = lim ( Da A(a)K 2N (a, a; t). (6) 

V / N^oo J 

It has already been proven that both AT(t) and T(t) Af_1 are Hilbert-Schmidt for all t > 
and thus AT(t/N) N is trace-class and its trace can be expressed in terms of the kernels of 
these operators. It is then a matter of inspection to conclude 

Tr (AT(t/N) N ^J = jDaA(a)K N (a,a;t) (7) 

and the theorem follows from equations (^) and (|7|). □ 

If liniTv K^(a, a; t) exists for almost all a, then an application of the dominated conver- 
gence theorem in ([]) shows that Tr^Ae^^j = lim^r fDaA(a)Kj^(a, a; t). Therefore, also 

odd lattice sizes can be included in the limit, for instance, if V is continuous (then is 
continuous and thus finite on the diagonal). 

In the second theorem we will work with temperatures having an imaginary part. The 
following proof relies on a simple property of the Hilbert-Schmidt operators given as Lemma 
IV. 1 in Section |V} To make the theorem easier to use, we will now state all the assumptions 



made on the potential explicitly. 

Theorem III. 5 Assume that A is an essentially bounded measurable function and V G L\ oc 
is a potential bounded from below for which ([Jj holds and 

jDae- tv{a) < oo, for all t > 0. 

Then for all complex (3 with Re/? > ; Ae~^ H is trace-class and 

Ti(Ae- pS ) = lim \DaA(a)K 2N (a,a](3). 

Proof: Let f N (J3) = fDa. A(a)K 2N (a, a; (5) and g{f3) = Tr^e"^). Both are then well- 
defined on the right half-plane: f^ is an 2iV-fold integral over an absolutely integrable 
function and, by the Schwarz inequality, self-adjointness of H and Theorem |111.4| , 

\g{(3)\ 2 < Tr(|l| 2 e- Re ^) Tr( e - Re/3 ^) < oo. 

In addition, g and all /jy are analytic on the right half-plane, as can be seen e.g. by a 
standard application of Morera's theorem — note that 

n 

where {E n } are the eigenvalues of H and {4>n} is a corresponding orthonormal basis. 



5 



Since Theorem [III.4| states that limAr^™ /jv(/3) = g((3) for positive (3, the theorem will 
follow from the Vitali convergence theorems if we can only prove that the sequence (/n) is 
uniformly bounded on every compact subset of the half-plane. 

An argument similar to one in the proof of Theorem [III.4| shows that Jn{P) — 

Tr^Af(^) 2N y Therefore, Lemma W^i establishes the bound \f N (j3)\ < \\A\\Tr (f *f ) N , 
where f = f (^). However, as f*f - ~~ 



H 0i 







Tr (T*T) N = Tr[T(Re/3/N)" ) < F (Ref3) < oo 



e 2jv v g jv u e 2jv 



N 



V 



by Corollary [111.2| and Theorem |111.4| . If S is a compact set in the right half-plane, it follows 
that Fo(ReS') is bounded since F Q was continuous. This proves that (/at) is uniformly 
bounded on every compact subset of the right half-plane and hence completes the proof. □ 

If V is continuous, then the limit can be taken by using all values of N: Define h^{(5) = 
JDslA(sl)Kn(sl, a; 0). By the remark after Theorem 111.4 , liniTv = g on the positive real 
axis. Also, 



\h 



2N+1 



Tr( AT( — - — ' 

V V 2N + 1' 



2N+1 



< \\A\\\\T\\Tt (T*T) N < 



F (Rep 2N ), 
uv 1 2N+l h 



and thus the sequence (h^) is uniformly bounded on the right half-plane. 



IV. PROOFS OF THE LEMMAS 



Proof of Lemma \III. 1\ : By the relation between geometric and arithmetic means, 
exp(— jj Yl!k=\ ^( x fc)) — lEf=i ex P( - ^( x *;))- Using this and the semigroup property 
of the integral kernels P, we arrive at the inequality 

A0v(a, b; t) < -e-^( b )p(a, b; t) + / Dc e~ tv ^- £ P(a, c; ^)P(c, b; (N -k)-). 



The terms in the last sum are explicitly 

n 

exp 



1 



N 



1 / N , 



2nt ^k{N - k) 

= P(a,c;i)P(c,b;i) 
N 



k=l 



N 



,a-cr + 
2t V fc 1 N-k 



Ic-bl 



N 



n-l 

N-k\ 2 



N-k 



^k(N -k)\ k 

n ~ l r fc |c-b| 2 

exp 



exp 



A- jfe |a- c| 
k 2~T 



A - /c 2t 



P(a,c;t)P(c,b;t) 



N 



n-l 
k \ 2 



X 



N 



y/k{N-k) \N-k 

n ~ l r A^-A;|a-c| 2 
exp 



exp 



|c — b| 

N — k 2~T 



2f 



(8) 



(9) 
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Assume now that n > 1 and consider first the case 1 < k < Since, t s e tr < (—) s for all 

— — 2 ' — \reJ 

positive r, s and t, the expression @ has an upper bound 



P(a,c;t)P(c,b;t) 



n-l 



y/k{N-k) 



e a — c 



On the other hand, if y < & < Af — 1, then the same reasoning can be applied to (||) leading 
to the same bound apart from the change of |a — c| to |c — b|. Since in both cases these 
upper bounds are positive, the terms under inspection are for all 1 < k < N — 1 bounded 
by 

n-l 

P(a, c; t)P(c, b; f) -L— -i f -— — T + * 



Now Xlfc^i 1 [77(1 — 77)] 2 is a lower Riemann sum for the integral ds— r^==y = 7r 
and, therefore for all iV, 

1 ^ JV 
Nf-[y/k(N-k) ~ 



n—1 n— 1 



and the Lemma follows if we choose C„ = 7r2 ri ~ 1 (n — l)^ _ e" t - . 

The case n = 1 is even simpler to prove, as then both of the exponentials in (|8]) can be 
majorized by one and the Lemma, with the choice C\ = 7T, follows from fllQ|) . □ 

Proofs of Corollaries \HI.Q and \III. ^ : Define the following two functions on (0, 00) 

I s (t):= y"pa|a| s e-* y(a) , ■/.(*):= J Dm |a| s P(0, a; t), 

where s is a real parameter. J s is finite for s > —n and, in fact, it is then continuous, since 
a change of variables and relying on n-dimensional spherical coordinates shows that 

ut) = m s/2 -^-. (ii) 

Under the assumptions of the theorem, also Iq and I r are everywhere finite. They are then 
also continuous, since for every converging sequence (tk) in (0, 00) dominated convergence 
can be used to prove linn^oo I s (tk) = J s (limt fc ). 

The rest is just a matter of bookkeeping when integrating both sides of the inequality 
([]). Using Jensen's inequality for r > 1 and the obvious relation (1 + |a|) r — |a| r < 1 for 
< r < 1 it is easy to prove that for r > and a, b G 1" 

|a| r < 2 r (|a-b| r + |b| r ). 

Applying this result, equation ( [LTD an d a little algebra shows that 

y'pa|a| r ir J v(a,a;t) < C {1 H^ I r (t) + C {2) t-^ I (t), 
jDaDb\a\ r K N (a,b;t) < C (3) I r (t) + C (4) ^I (t), 
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where the constants depend only on n and r. The functions on the right hand sides are 
continuous and can therefore be chosen as F r and G r . 

Since e~ tv ^P(c, b; t) < (2nt)~ n / 2 , Corollary [III.3| is an immediate consequence of (|Tl~D 
and Lemma [111.1| . □ 

Lemma IV. 1 If A is bounded, T is Hilbert-Schmidt and N > 1, then 



Tr [AT 



2N 



< \\A\\Tt(T*T) 



N 



Proof: By the Schwarz inequality, 



Tr ( AT 



2N 



< Tr 



(AT N )*(AT 



Tr 



^rpN\*fj-\N 



and since Tr 



(AT N )*(AT 



Tr 



< 



! Tr 



we get 



Tr [AT 



2N 



< WAWTt 



Using a method similar to the so called min-max principle, K. Fan has shown0 that 
for any finite-dimensional operator (i.e. matrix) B, Tr (B N )*B N < Tr (B*B) N . The finite- 
dimensional operators are dense in T 2 = the Hilbert space of Hilbert-Schmidt operators with 
the inner product (■, -)hs- Since the product, the inner product and the norm are continuous 

on 1~2 and Tr 



(B N )*B N 



B N \\l s , Tr (B*B) N = (B*B,(B*B) N - 1 ) m , it is obvious that 



the result of Fan is valid for every Hilbert-Schmidt operator. 



□ 



V. DISCUSSION 



The main results of this paper are the two theorems, HL4| and |IL5[ The first one justifies 



the faith of the lattice community that also expectation values of polynomials are measurable 
with the lattice Monte Carlo methods if the potential increases sufficiently fast at infinity. 
The second one establishes the possibility of using the lattice integrals for computations of 
the expectation values even when the temperature has an imaginary part. 

There are a few shortcomings of the present results. First, the exact class of potentials 
for which the integral kernels can be used was not determined. This is not a serious 
problem since, by accepting the necessity for a double limit, statistics given by any potential 
in L\ oc can be approached with the lattice integrals by replacing V by a sequence 14 G C°° 
which converges to it pointwise almost everywhere. The second shortcoming is the class of 
observables in Theorem |111.5| : the present proof cannot be easily extended to include also 
polynomial observables for polynomial potentials, although it is very natural to expect the 
integrals to converge to the right limit by Theorem [111.4| . Finally, the necessity of using only 
even number of lattice points has no natural explanation and as such it is expected to be an 
artifact of the present approach. Indeed, as has been noted earlier, the requirement of even 
lattice sizes can be dropped for continuous potentials. 
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